Seasonal dynamics of ammonia-oxidizing bacteria but not archaea influence soil nitrogen cycling in a semi-arid agricultural soil

Nitrification, a key pathway of nitrogen (N) loss from agricultural soils, is performed by ammonia-oxidizing bacteria (AOB) and archaea (AOA). We examined the seasonal dynamics (2 years) of ammonia oxidizer gene abundances across a gradient of soil carbon (C) and N in a semi-arid soil after 8 years of tillage and crop residue treatments. AOB was more dominant than AOA in the surface soil, as AOA were undetected in 96% of samples. Seasonal variation in AOB abundance was related to substrate availability; AOB gene copy numbers increased at the end of the growing season (during summer fallow) following higher concentrations in dissolved organic matter soil water. This suggests increased co-location between AOB and substrate resources in pores still filled with water as the soils dried. AOB was however not statistically related to soil ammonium concentrations, soil water content, rainfall or temperature. Organic matter inputs enhanced AOB abundance independent of seasonal variation. AOB abundance was greatest in autumn and immediately preceding the start of the growing season, and coincided with elevated soil nitrate concentrations. The growth of the AOB population is likely to contribute to increased risk of N loss through leaching and/or denitrification at the start of the crop growing season following summer fallow.

www.nature.com/scientificreports/ with a meta-analysis of agricultural soils demonstrating N fertilizer form and soil pH affecting the response of AOA to N fertilization 20 .
Factors that regulate temporal variation in populations are often different from factors that regulate spatial fluctuations of populations 22 . Temporal variation in ammonia oxidizer gene abundance may be regulated by season (e.g. temperature and rainfall), soil organic matter (OM), and changes in NH 4 + availability due to N fertilizer additions and periods of soil OM mineralization 8,[23][24][25] . Taylor et al. 24 speculated that increases of AOA amoA abundance (but not AOB) in spring of a fallow Oregon agricultural soil (with a dry-summer subtropical climate) was due to their utilization of another energy-generating metabolism besides ammonia oxidation (i.e. heterotrophy or mixotrophy). These authors also suggested that increased AOB amoA gene abundance in cropped soil at the same site was in response to increased NH 4 + availability from N fertilizer additions 24 . In a semi-arid soil in Israel, greater AOA abundance than AOB in summer was hypothesized to be due to low NH 4 + availability, as a consequence of low soil water content causing low mineralization rates 26 . Ammonia-oxidizing bacteria in the same soil may have also been limited by elevated temperature to a greater degree than AOA: in winter AOB shifted to greater dominance than AOA in the same soil 26 . However in agricultural soils, seasonal changes in the relative abundance of AOA to AOB may also be attributed to changes in nutrient availability in response to crop growth 8 .
The effect of soil OM content on amoA gene abundance of ammonia oxidizing bacteria and archaea varies. Many archaeal ammonia oxidizers appear to be able to take up organic C as an energy source (show heterotrophic and mixotrophic growth), so would be expected to be more abundant in soils with greater OM contents 27 . Bacterial ammonia oxidizers are limited to autotrophic growth, but increased supply of ammonia substrate from enhanced mineralization in higher OM soil might also increase AOB population size. However, studies have observed that increased OM or increased soil C and N pools have, for example, a greater AOA:AOB ratio 28 ; no effect on AOB gene abundance 29,30 ; or are negatively correlated to AOA gene abundance 30 . Temporal effects can further interact with spatial effects, such as an observed doubling of AOB gene abundance at the start of the growing season due to organic inputs that increased soil C and N pools, but not at other times of the year 31 ; and a negative correlation between AOB but not AOA abundance and soil OM content changes with season 26 . These temporal changes may be due to the changes in soil water content and subsequent OM mineralization and NH 4 + availability 26 . Drylands, including semi-arid and arid regions, are important for agriculture, comprising one third of the world's agriculturally productive lands and supporting nearly half of the global population 32,33 . Semi-arid agricultural soils in the grain-growing region of south-western Australia are generally acidic, sandy, low in soil OM content, and are annually exposed to soil temperatures greater than 40 °C during summer [34][35][36] . Rain in this region mainly falls during winter (the growing season), while agricultural soil is often left fallow during the dry summer (Fig. 1a). Soil NH 4 + concentrations are usually less than 5 µg N g −1 dry soil except for short periods when NH 4 + -based fertilizers are applied to agricultural soil during the winter growing season 35,36 . Archaeal ammonia oxidizers would be expected to have greater abundance than AOB in these soils, due to their better resistance to elevated temperatures, low pH and low NH 4 + availability 18,26,37 . Studies of surface soils in this region however have shown that AOB are similar to or more abundant than AOA 14,17,38,39 . As these studies sampled soil at one point in time, it is not yet clear whether this dominance of AOB occurs across all seasons, and how AOB and AOA populations vary dynamically with time in response to soil and environmental conditions.
Here we aimed to understand how soil OM that had been built by agricultural management practices in one semi-arid soil affected amoA gene abundance, and whether these interacted with season to modify population sizes of ammonia oxidizers. The overall objective of this two-year field-based study was to assess the temporal population dynamics of ammonia oxidizers so as to better understand potential N loss mechanisms in semi-arid environments dominated by winter rainfall. Specifically, this study investigated: (i) seasonal variation in ammonia oxidizer abundance; (ii) which soil and environmental factors regulated the relative abundances of AOB and AOA; and (iii) whether increased soil OM modified seasonal variation in ammonia oxidizer abundance. Five field-based OM treatments were chosen that were designed to alter soil C and N availability: no tillage ("No Till"), no till with burnt stubble ("Burnt"), tillage ("Tilled"), tillage with OM additions ("Tilled + OM"), and tillage with OM additions that ceased four years before this study commenced ("Tilled + OM Run-Down"). We hypothesized that there would be a greater abundance of AOB than AOA, AOB abundance would vary temporally in response to dissolved OM, and increasing soil OM would increase AOA relative to AOB throughout the year.

Results
Environmental conditions. Daily rainfall over the study period ranged from 0 to 25.2 mm (Fig. 1a).
Growing season rainfall (April-October) was greater in 2011 (283 mm) than in 2010 (144 mm). Daily minimum soil temperature at 5 cm depth ranged from 6.6 °C in July 2011 to 30.4 °C in January 2011 (Fig. 1a). Daily maximum soil temperature at 5 cm depth ranged from 9.2 °C in August 2010 to 41.7 °C in January 2011 (Fig. 1a).
Soil water content at the time of soil collection followed a similar pattern to rainfall, increasing during the winter growing season (May-November) and at a minimum during summer (mean across OM treatments < 0.01 g g −1 dry soil; October 2010-May 2011; Fig. 1b). Mean maximum soil water content across OM treatments at sampling was greater in 2011 (0.09 g g −1 dry soil in August) than in 2010 (0.05 g g −1 dry soil in July-August; p < 0.0001). Tilled + OM soil had greater soil water content than the other OM treatments on two occasions in winter (July 2010, August 2011; p < 0.01), and greater soil water content than the No Till, Burnt and Tilled soils on one occasion in late spring (November 2011; p < 0.05; Fig. 1b).
Dissolved organic C and microbial biomass C. Dissolved organic matter (DOM) was analyzed as dissolved organic C (DOC). Mean DOC ranged from 60 to 217 µg C g −1 dry soil (Fig. 2a)  www.nature.com/scientificreports/ water (Fig. 3a). During the dry summer months mean DOC concentration in soil solution averaged 113 mg C g −1 water (December 2010) compared to during winter where it is diluted to below 26 mg C g −1 water (May-September; Fig. 3a). However, DOC per gram of dry soil was at a maximum in December 2010 (summer fallow) and at a minimum in August 2011 (winter, mid-growing season; Fig. 2a). Although DOC in soil solution was statistically different between some OM treatments on each sampling date, the range between these treatments was very small (maximum statistically significant difference was 1.1 mg C g −1 water between Tilled and Tilled + OM soils on 19th Sep 2011; p < 0.001) compared to the variation of DOC over time. Dissolved organic C in Tilled + OM soil was greater than No Till, Burnt and Tilled soils at all sampling times (p < 0.01), and was greater than Tilled + OM Run-Down from July-December in 2010 (mid-winter to mid-summer), and from September-December in 2011 (spring to mid-summer; p < 0.01; Fig. 2a). Dissolved organic C in No Till, Burnt and Tilled soils were not different from each other (p > 0.05).
Mean microbial biomass C (MBC) ranged from 46 to 489 µg C g −1 dry soil (Supp. Fig. 1a). Microbial biomass C followed similar patterns with time in all OM treatments, but the Tilled + OM soil had significantly greater MBC than the other OM treatments (p < 0.05), except in May 2010 (autumn) and June 2011 (winter).
Mean NH 4 + in soil solution peaked in summer when soil water content was lowest (December 2010; p < 0.0001; Fig. 3b). However, mean NH 4 + per gram of dry soil in No Till, Burnt and Tilled soils peaked in winter (July 2010 and 2011), coinciding with N fertilizer inputs (Fig. 2b). Similar to DOC, NH 4 + in soil solution was statistically different between some OM treatments on each sampling date, but these differences were small compared to seasonal variation in soil solution NH 4 + (maximum statistically significant difference was 0.7 mg N g −1 water between Burnt and Tilled + OM RD soils on 20 Jun 2011; p < 0.0001). Surprisingly, NH 4 + concentrations in the Tilled + OM soil did not respond to N fertilizer inputs, however NH 4 + concentrations in Tilled + OM soil were greater than the other OM treatments during summer (October 2010-May 2011; p < 0.05).
Nitrate in soil solution was on average 816% greater in December 2010 (summer) and May 2011 (autumn) than at other times of the year (p < 0.05; Fig. 3c). Nitrate concentration in soil solution decreased significantly at the beginning of autumn growing season rains (May-June 2011; p < 0001 for all OM treatments). Soil www.nature.com/scientificreports/ NO 3 − concentrations in all OM treatments also increased from the time of harvest by up to 1069% (p < 0.0001), in response to a series of summer rainfall events (compare Figs. 1a, 3b). Soil NO 3 − peaked at the start of each growing season, and then declined to a minimum by the end of winter (around September; p < 0.0001; Fig. 2c). Mean NO 3 − concentrations in Tilled + OM soil and Tilled + OM Run Down were greater than the other OM treatments in soil (p < 0.001) and in soil water over summer/autumn (May and December; p < 0.05).
Mean potentially mineralizable N (PMN) ranged from 10 to 98 µg N g −1 dry soil (Supp. Ammonia-oxidizer gene abundance. Archaeal amoA gene abundance was below detection limits for 96% of samples (144 out of 150 samples; data not shown). Of the six samples in which AOA were detected, five were from No Till soil at different time points across all seasons sampled, and one was from the Tilled + OM Bacterial amoA gene abundance ranged from 2.35 × 10 7 to 1.11 × 10 9 gene copies g −1 dry soil (Fig. 4), and was detected in all samples. AOB varied seasonally in a similar way in each OM treatment. Bacterial amoA gene abundance generally decreased from the middle to the end of each growing season and was at a minimum in spring to early summer. Bacterial amoA gene abundance increased on average by 839% over summer to a maximum at the start of the growing season (from Dec 2010 to May 2011; p < 0.0001). When averaged across seasons, Tilled + OM soil had greater bacterial amoA gene abundance (mean of 4.91 × 10 8 gene copies g −1 dry soil) than the other OM treatments (p < 0.05), and Tilled + OM Run-Down had greater bacterial amoA gene abundance (mean of 3.63 × 10 8 gene copies g −1 dry soil) than Burnt and Tilled soils (means of 2.15 × 10 8 and 2.17 × 10 8 gene copies g −1 dry soil respectively; p < 0.05), but was not different from No Till soil (mean of 3.07 × 10 8 gene copies g −1 dry soil; p > 0.05).

Relationships between bacterial amoA gene abundance and other variables. Separate linear
regressions of each variable were better related to bacterial amoA gene abundance (had the lowest Akaike Information Criteria) than multiple linear regressions using more than one variable. Bacterial amoA gene abundance   Fig. 2), but no relationship with cumulative 30-day rainfall, mean daily maximum or minimum soil temperature in the 30 days prior to sampling, soil water content at collection, NH 4 + (in soil or in soil water), NO 3 − in soil water, DOC in soil water or PMN (p > 0.05). The potential relationship between AOB and each significantly correlated variable was calculated by multiplying the coefficient by the maximum value measured for that variable. Variables with linear regression coefficients in order of greatest to least correlation to AOB were DOC, MBC, total inorganic N and NO 3 − . Principal component analysis (PCA) showed that samples from the same collection date generally grouped across the biplot of principal components 1 (PC1) and 2 (PC2), roughly in line with PC1 as rainfall, soil water content and soil temperature varied (and consequently DOC, NH 4 + and NO 3 − in soil water varied; Fig. 5, collection dates shown by symbol shapes). Samples from the Tilled + OM soil, and Tilled + OM Run-Down to a lesser extent, had greater scatter across the biplot in line with PC2, compared to the No Till, Burnt and Tilled soils, which were grouped towards positive scores of PC2 (Fig. 5, OM treatments shown by symbol colours).
Five principal components were required to explain greater than 90% of the variance in soil biochemical and environmental properties (Supp. Table 1). PC1 explained 50.2% of total variance in the data (Supp. Table 1), and separated most strongly cumulative 30-day mean minimum and maximum soil temperatures (positive loadings)  www.nature.com/scientificreports/ from cumulative 30-day rainfall and soil water content at collection (negative loadings; Fig. 5; Supp. Table 2). Otherwise, PC1 separated DOC in soil water, NH 4 + in soil water and NO 3 − in soil water (positive loadings near minimum and maximum temperature) from MBC, PMN and logged AOB gene abundance (with slightly negative loadings). PC2 explained 15.9% of total variance in the data (Supp. Table 1), and separated cumulative 30-day rainfall, soil water content at collection, DOC in soil water and NH 4 + in soil water (negative loadings) from MBC, logged AOB gene abundance, and PMN (most strongly negative loadings; Fig. 5; Supp. Table 2). Otherwise, cumulative 30-day mean minimum and maximum temperatures and NO 3 − in soil water were grouped together on PC2 (slightly positive loadings).

Discussion
Ammonia-oxidizing bacteria rather than AOA appear to regulate nitrification processes in the surface layer of this semi-arid agricultural soil across all seasons. This agrees with previous studies from this region 14,38,39 , but is in contrast to other semi-arid environments and agricultural soils, which have found that AOA are often dominant to AOB, particularly in soils with low pH 23,[39][40][41] . Here, we hypothesize that AOB were present but AOA were not detectable in the surface soil (96% of samples were undetected) as a consequence of the differing ecophysiologies and trace mineral requirements of each microorganism. Electron transfer enzymes in AOA are based on copper, in contrast to the iron-dependent ammonia oxidizing systems of AOB 13 . Soils in this Western Australian semi-arid region are ancient and highly weathered, so are rich in iron but deficient in copper 42,43 . Jenkins et al. 14 observed a significant relationship between AOA and copper content in agricultural soils across Australia, suggesting that copper limitation of AOA is an important factor regulating niche differentiation between ammonia oxidizing organisms in soils in this region. Consequently, AOB is a target group for decreasing N losses, via leaching and denitrification, and improving N use efficiency in semi-arid regions containing highly weathered soils.
Seasonal variation in AOB abundance was primarily related to substrate concentration in pore water. Dissolved organic matter (here measured as DOC) and NH 4 + in soil solution showed clear peaks in late winter at the end of the crop growing season (Fig. 3a,b), before an increase in AOB abundance during summer and autumn fallow and subsequent production of NO 3 − (Figs. 3c, 4). Mineralization is inferred to be stimulated by crop plants (the source of DOM), particularly at the end of the growing season. Active nitrification and rapid depletion of NH 4 + pools follows, as indicated by NO 3 − dominating inorganic N pools in this soil (Figs. 2, 3  www.nature.com/scientificreports/ soil solution volume in this semi-arid soil appear to have as much of an effect on substrate availability to microorganisms as the production of those substrates, through the concentrating effect of soil drying (i.e. physically locating the substrate into pores and water films where microorganisms are located). The poor relationship between bacterial amoA gene abundance and rainfall, soil temperature, and soil water content was unexpected. Other authors have observed that climate variables can be as important, if not more, than agricultural management effects on soil and microbial processes in semi-arid soils 44,45 . For example, Wang et al. 41 showed that AOB abundance was more sensitive to short-term N inputs and warming than AOA in a semi-arid grassland 26 . The lack of a relationship between AOB abundance and environmental conditions in the present study might be explained by a disconnect between longer-term population dynamics and shorter-term environmental events. Although directly influenced by environmental conditions such as soil water and temperature, abundances of functional genes are less dynamic than many environmental and edaphic properties, so functional gene abundances are likely to reflect the longer-term variation of those properties 46 . There was no linear correlation between AOB and climate characteristics, however the PCA showed that the variation in the AOB abundance, soil properties and climate data was partly explained by rainfall, soil water content and temperature, and clearly grouped by sampling date (Fig. 5, PC1 in Supp. Table 2). These findings suggest that seasonal dynamics of AOB populations are not as yet well understood, but are due to many interacting factors that will consequently influence risk of N loss through the production of NO 3 − . Organic matter inputs to soil have positively but indirectly influenced bacterial amoA gene abundance (irrespective of the time of year) due to the enhanced supply of either C or N substrates. The effect of soil OM content on AOB in the present study did not support our hypothesis and varied from observations in other semi-arid soils. For example, whilst bacterial amoA gene abundance was positively correlated with soil OM in our study, it was negatively correlated with soil OM in summer and winter over six years in a semi-arid shrubland soil 26 . This negative correlation would be expected to occur if soil OM stimulates heterotrophic microorganisms to compete more successfully for NH 4 + , decreasing N substrate availability and thus growth of autotrophic ammonia oxidizers. The positive effect of increased soil OM pools on AOB abundance in the present study could be attributed to the associated increases in N supply through DOM to autotrophic ammonia oxidizers. Previous work on the same soil showed that OM inputs increased gross N mineralization, immobilization as well as nitrification rates 47 , suggesting that increased soil OM stimulates both heterotrophic and autotrophic microbial populations and activity.
Bacterial ammonia oxidizer growth during summer fallow is likely an important factor determining increases in soil NO 3 − pools and consequently risk of N loss in this semi-arid soil. Soil NO 3 − pools are a balance of inputs and outputs. Nitrate inputs into the surface soil (0-10 cm) appear to continue as summer progresses with maintenance of NO 3 − concentrations in soil solution (Fig. 3c) and increased total soil NO 3 − (Fig. 2c). Nitrate inputs over the dry season has also been observed in the same and similar semi-arid soils in response to summer rainfall events that stimulate N mineralization 35,36,48 . In contrast to inputs, outputs of NO 3 − in this region are likely to be low during summer fallow, with plant uptake minimal or non-existent due to the absence of plants, and it is not until the start of the growing season rains that any significant deep drainage and NO 3 − leaching usually occurs 49 . A significant decrease in concentration of NO 3 − in soil solution was observed in the present study in early autumn (Fig. 3c), likely due to leaching and dilution by rainwater. Correlation of bacterial amoA gene abundance and NO 3 − pool size in soil suggests that increases in NO 3 − during summer are not due merely to N mineralization flushes and lack of inorganic N uptake, but also due to growth of the AOB population, strongly coupled to activity of nitrite-oxidizing bacteria. Bacterial amoA gene abundance was less in summer than in winter in another semi-arid soil from Israel, which was attributed to the sensitivity of AOB to high temperatures, at that site ranging from 23 to 35 °C in summer 26 . By contrast, bacterial amoA gene abundance in the present study continued to increase over summer, even though soil temperatures exceeded 40 °C (Fig. 1a). This may be because AOB were sampled more often, so finer scale population dynamics could be observed in this soil. An alternative explanation is that AOB in the present study region have adapted to periods of high temperature and low water availability, allowing them to survive then continue to grow when conditions become favourable, even if only briefly due to transient water pulses from summer rainfall events.
Predicted changes in climate for this region are likely to exacerbate the effects of AOB population dynamics and NO 3 − production, particularly during summer fallow. Since the 1970s, mean annual rainfall in Southern Hemisphere semi-arid regions has been decreasing, especially during autumn and winter at the start of the growing season 50,51 . Summer rainfall events however are increasing 52 , which will likely enhance inorganic N supply through OM decomposition when there is no plant N uptake from fallow soil 48,53 . This inorganic N is then at risk of loss if it is in excess of microbial demand and is subjected to nitrification. Managing undesirable losses of N from soil in response to summer rainfall events change in this semi-arid region is likely to become increasingly challenging, unless the activity and growth of the AOB population can be controlled (or excess inorganic N utilized by plants or soil microorganisms before it is lost from the soil).
In conclusion, seasonal variation of bacterial amoA gene abundance in this soil is contributing to NO 3 − production and therefore the risk of N loss, especially as summer fallow progresses. Archaea however are unlikely to be influential drivers of ammonia oxidation as archaeal amoA gene abundances were predominantly below detection limits in the surface layers, possibly due to copper limitation of archaeal ammonia oxidation enzyme systems. Increased soil OM levels by additional crop residue inputs positively influenced bacterial amoA gene abundances, but did not modify seasonal variation in ammonia oxidizer abundance. Increases in AOB abundance over summer fallow followed a pulse of DOM and NH 4 + in soil solution at the end of the winter growing season. As expected, bacterial amoA gene abundance was not statistically related to NH 4 + concentration in soil, but surprisingly was also not related to soil water content, 30-day rainfall or mean soil temperature, suggesting that ammonia oxidizer populations in these semi-arid soils are regulated by longer-term changes in climate and substrate supply. OM was applied after each lupin crop, and before seeding of wheat. Seven years after the study commencement, soil organic carbon (SOC) contents ranged from 0.91 to 1.38%, and total N contents ranged from 0.09 to 0.13% (Table 2). Two crops of winter wheat were grown during the present study, and fertilized at rates depending on expected growing season rainfall and projected potential yields (1. Soil was collected (0-10 cm depth) on ten occasions from each plot between May 2010 and November 2011. Within each plot, 25 cores (5.3 cm diameter × 10 cm depth) were sampled using a zigzag pattern. Subsamples of soil from each plot were frozen for DNA analysis of AOA and AOB on return from the field. Remaining soil was sieved to < 2 mm and stored field-moist at 4 °C. Soils were analysed for gravimetric soil water content and inorganic N at field water content. All soils were adjusted to 45 to 55% water holding capacity 56 , and then preincubated for 7 days at 25 °C before analysis for MBC, DOC and PMN. On the other sampling dates there was no need to first wet-up the soil.
Microbial biomass C and dissolved organic matter. Microbial biomass C was measured using the chloroform fumigation-extraction method 56 . Fresh soil (10 g) was fumigated with chloroform (stabilized by amylene) under vacuum in the dark for 24 h at 25 °C. Fumigated samples, with replicate non-fumigated samples, were extracted by shaking for 1 h with 40 mL of 0.5 M K 2 SO 4 . Filtered extracts (Whatman No. 42) were frozen until further analysis. Fumigated and non-fumigated extracts were analysed using an OI Analytical Aurora Table 2. Properties of field organic matter treatments (0-10 cm depth) at start of present study, 7 years after treatments were imposed. Values are a mean of three field replicates ± standard error of the mean. Organic matter treatments with the same letter are not significantly different (p > 0.05). OM organic matter. # Bulk density determined using the intact core method with 3 cores of 7.35 cm diameter by 10 cm depth 66 . § pH was determined on air-dry soil in 0.01 M CaCl 2 with a 1:5 soil:extract ratio, after shaking for 1 h, and while stirring the soil suspension 67 . ‡ Electrical conductivity was determined on air-dry soil in water with a 1:5 soil:water ratio 67 . † Total carbon and nitrogen were determined by high-temperature combustion of finely ground air-dry soil using an Elementar Vario MACRO CNS elemental analyzer (Hanau, Germany) 67 . ∞ Trace elements copper and iron were determined by atomic absorption spectroscopy after extraction with diethylene-triamine-pentaacetic acid (DTPA) solution (ratio of 1:2) for 2 h 67 . . Both NH 4 + and NO 3 − were expressed per gram of oven-dry soil and per gram of soil water as described above.
Potentially mineralizable N was determined by anaerobic incubation 60,61 . Fresh soil (20 g) was incubated in 80 mL of deionized water for 7 days at 40 °C. K 2 SO 4 was added (6.97 g) to adjust the soil solution to 0.5 M K 2 SO 4 , then samples were shaken for 1 h. Replicate 20 g samples of non-incubated soil in 80 mL of 0.5 M K 2 SO 4 were also shaken. Incubated and non-incubated samples were filtered (Whatman No. 42) then frozen until further analysis. Samples were analysed for NH 4 + on an autoanalyser, as described above.
Nucleic acid extraction and qPCR. In March 2014, DNA was extracted from 700 mg subsamples of field-moist soil 62  Archaeal amoA primers were Arch-amoAF (5′-STAA TGG TCT GGC TTA GAC G-3′) and Arch-amoAR (5′-GCG GCC ATC CAT CTG TAT GT-3′) amplifying a fragment length of 635 bp 63 . Archaeal amoA gene abundance was quantified using an Applied Biosystems ViiA™ 7 (Carlsbad, CA, USA). Each 10 µL qPCR reaction contained 5 µL of Power SYBR Green PCR Master Mix (Applied Biosystems, Warrington, UK), 0.1 µL each of the specific forward and reverse primers at 10 µM, 1 µL of BSA, 1 µL of template DNA and 2.8 µL of water. Thermocycling conditions were: 94 °C for 10 min; then 40 cycles of 94 °C for 1 min, 52 °C for 1 min, 72 °C for 1 min and 78 °C for 1 min; followed by a melt curve to establish a single product amplification. Fluorescence data was collected at the 78 °C stage.
Standard templates used to determine gene copy numbers in the qPCR reactions were cloned plasmids, as described in Barton et al. 64 . Samples were tested over a series of dilutions to determine if there was inhibition, and further analysis was completed using the dilution that produced the highest copy number. Standard curves generated in each reaction were linear over four orders of magnitude for AOB (10 3 -10 6 gene copies) and six orders of magnitude for AOA (10 3 -10 8 ) with r 2 values greater than 0.98. Amplification efficiencies ranged from 79 to 8% for AOA amoA quantification and 94-111% for quantification of AOB amoA. Statistical analysis. Statistical differences between OM treatments for the basic soil properties were determined using ANOVA with associated TukeyHSD post hoc tests in R version 3.1.0 (R Foundation for Statistical Computing, Vienna, Austria). Bacterial amoA gene abundance was initially log transformed before all statistical analyses. Statistical significances of OM treatment with time for each soil property were evaluated by mixed linear models estimated by restricted maximum likelihood (REML), using PROC MIXED in SAS version 9.3 (SAS Institute Inc., Cary, NC, USA), because replication was at the highest level (i.e. replicates were field plots). Random effects of field block and field block interacting with OM treatment were included in the model. For this analysis, NH 4 + in soil and in soil water, PMN, MBC in soil and soil water and DOC in soil and soil water were also log transformed, while total inorganic N, NO 3 − in soil and soil water and AOB abundance were square root transformed to ensure a normal distribution. When there was a significant interaction between OM treatment and date, and multiple pairwise comparisons were needed to support results statements, the reported P value is the maximum P value of the relevant pairwise comparisons.
Backwards multiple linear regression was also conducted using PROC MIXED. Random effects of field block and field block interacting with OM treatment were included in the model. The best regression model was evaluated by minimizing the Akaike Information Criteria 65 . PCA was carried out using prcomp in R version 3.1.0 with scaled and centred data.
Experimental protocols. All plant experiments were in compliance with relevant institutional, national, and international guidelines and legislation.

Data availability
The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.